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The geometric scaling property observed in the HERA data at small x, that the deep inelastic 
scattering (DIS) total cross-section is a function of only the variable Q 2 x x , has been argued to 
be a manifestation of the saturation regime of QCD and of the saturation scale Q 2 a (x) ~a; _A . We 
show that this implies a similar scaling in the context of diffractive DIS and we observe, for several 
diffractive observables, that the experimental data from HERA confirm the expectations of this 
scaling. 
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I. INTRODUCTION 

o : 

Deep inelastic scattering (DIS) is a process in which a virtual photon is used as a hard probe to resolve the small 
distances inside a proton and study its partonic constituents: quarks and gluons that obey the laws of perturbative 
i— I \ QCD. When probing with a fixed photon virtuality Q 2 3> Aq C£) , and increasing the energy of the photon-proton 
I ^ ■ collision W, the parton densities inside the proton grow. Eventually, at some energy much bigger than the hard scale, 
corresponding to a small value of x = Q 2 /W 2 , one enters the saturation regime of QCD the gluon density becomes 
so large that non-linear effects like gluon recombination become important, taming the growth of the parton densities. 
The transition to the saturation regime is characterized by the so-called saturation momentum Q s (x) — Qo x~ A//2 . 
, This is an intrinsic scale of the high-energy proton which increases as x decreases. Qo ~ Aqcd, but as the energy 
increases, Q s becomes a hard scale, and the transition to saturation occurs when Q s becomes comparable to Q. The 
^\ , higher Q 2 is, the smaller x should be to enter the saturation regime. 

Although the saturation regime is only reached when Q s ~ Q, observables are sensitive to the saturation scale 
already during the approach to saturation [2] when Aqcd^-Qs^Q- For inclusive events in deep inelastic scattering, 
this feature manifests itself via the so-called geometric scaling property: instead of being a function of Q 2 /Qq and 
x separately, the total cross-section is only a function of r = Q 2 JQ 2 {x), up to large values of r. Experimental 
, measurements of inclusive DIS are compatible with that prediction [3j . 

Part of the DIS events are diffractive, meaning that the proton remains intact after the collision and that there 
is a rapidity gap between that proton and the rest of the final-state particles. Such events are expected to be 
much more sensitive to the saturation regime of QCD than the inclusive ones Q and our goal in this letter is to 
extend the geometric scaling property to diffractive DIS and to compare the resulting predictions with the available 
experimental data. We shall consider several diffractive observables, focusing on inclusive hard diffraction and vector- 
meson production. 

In the saturation regime of QCD, contributions to the cross-sections growing like Q s /Q are important. The 
leading- twist approximation of perturbative QCD, in which Q 2 is taken as the bigger scale, cannot account for 
such contributions, and therefore is not appropriate to describe the small-x limit of deep inelastic scattering. As 
leading-twist gluon distributions cannot be used to compute cross-sections, the dipole picture of DIS [j| has been 
developed to describe the high-energy limit. It expresses the hadronic scattering of the virtual photon through its 
fluctuation into a color singlet qq pair (or dipole) of a transverse size r~l/Q. The dipole is then the hard probe that 
resolves the small distances inside the proton. 

The dipole picture naturally incorporates the description of both inclusive and diffractive events into a common 
theoretical framework 0, , as the same dipole scattering amplitudes enter in the formulation of the inclusive and 
diffractive cross-sections. This will be recalled in Section II and will allow us to extend the geometric scaling property 
to diffractive observables, as detailed in Section III. Finally, Section IV is devoted to comparison with experimental 
data and Section V concludes. 
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II. THE DIPOLE PICTURE OF DEEP INELASTIC SCATTERING 



We focus on diffractive DIS: 7*75— s-Xp. The proton gets out of the 7* — p collision intact, and there is a rapidity gap 
between that proton and the final state X whose invariant mass we denote Mj. We recall that the photon virtuality 
is denoted Q 2 , and the 7*— p total energy W. It is convenient to introduce the following variables: 

r>2 Q2 

I 3 = ^2 I n,r2 ' X P = X /l 3 ■ C 1 ) 



Q 2 + W 2 ' Q 2 + M\ 

The 7*— p total cross-section a1 f^' X is usually expressed as a function of x and Q 2 , while the diffractive cross-section 
da~ 1 di p J* Xp I d(3 is expressed as a function of j3, xp, and Q 2 . Note that the size of the rapidity gap in the final state is 
ln(l/a;p). 

To compute those cross-sections in the high-energy limit, it is convenient to view the process in a particular frame 
called the dipolc frame. In this frame, the virtual photon undergoes the hadronic interaction via a fluctuation into 
a colorless qq pair, called dipole, which then interacts with the target proton. The wavefunctions if)^ ^ {z,v;Q 2 ) 
describing the splitting of a virtual photon with polarization A into a dipole are well known. The indices a and (3 
denote the spins of the quark and the antiquark composing the dipole of flavor /. The wavefunctions depend on Q 2 , 
the fraction z of longitudinal momentum (with respect to the "f*—p collision axis) carried by the quark, and the 
two-dimensional vector r whose modulus is the transverse size of the dipole (transverse coordinates are obtained from 
a Fourier transform of transverse momenta). Formulae giving the funcions ip\' a ^ can be found in the literature (see 
for instance 0). In what follows, we will need the functions which describe the overlap between two wavefunctions 
for splitting into dipoles of different transverse size r and r' : 



b{(z, r, r'; Q 2 ) = N c £ [</>{' Q/3 (z, r'; Q 2 )] * ^{z, r; Q 2 ) . (2) 

a/3 



For a transversely (T) or longitudinally (L) polarized photon, these functions are given by 



{{z, r, r'; Q 2 ) = ~^e) [{z 2 + (1 - z) 2 )e) — Kl {s f \v\)Ki{e M') + m}K (s f \r\)K (e f \r\') 



(3) 



^{z^v<-Q 2 ) = ^^e 2 AQ 2 z 2 {l~z) 2 K {eM)K»(eM')- (4) 

In the above, e/ and mj denote the charge and mass of the quark with flavor / and e 2 =z(l — z)Q 2 +m'j. 

The total cross-section a] a f^ X . Via the optical theorem, it is related to the elastic scattering of the virtual 
photon. In the dipole frame, this scattering happens as follows: at a given impact parameter b, the photon splits into 
a dipole with a given size r which scatters elastically off the proton and recombines back into the photon. Therefore 
the overlap function $^ 7 which enters in the computation of the total cross-section is 

$fMz,|r|;Q 2 )=^{(z,r,r;Q 2 ) . (5) 
/ 

The total cross-section is then given by (for fixed impact parameter b): 
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,Q 2 )-2 f d 2 r [ dz J2 *r 7 *(^|r|;Q 2 ) T g(? (r,b;x) 



<jlh {x,Q*) = 2 J d'r dz ^ $1 (:. r :Q 2 ) r„„(r.b:.r) ((>) 

where the function T q q(r,h; x) is the elastic scattering amplitude of the dipole of size r off the proton at impact 
parameter b. It contains the x dependence, reflecting the fact that in our frame, the proton carries all the energy and 
is therefore evolved up to the rapidity ln(l/x). In the high-energy limit id we are considering here, T q q does not 
depend on z. 

The diffractive cross-section da 1 di jJ* Xp / d(3. In the amplitude, the photon splits into a dipole of size r which 
scatters off the proton and dissociates into a final state of invariant mass Mx- The same happens in the complex 
conjugate amplitude, except that the dipole size r' is different from r. Indeed, the final state is characterized by a 
particular value of Mx (or equivalently (3), corresponding to a particular momentum of the quark-antiquark pair. 



In coordinate space, this imposes two different dipole sizes in the amplitude and the complex conjugate amplitude, 
therefore the functions <fi x (z, r, r'; Q 2 ) (see (j2J) enter in the computation of the diffractive cross-section: 
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7*p — *Xp 
diff 

d 2 b dp 



-(/3,z P ,Q 2 ) 
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J d 2 r J d 2 r' J dzz{\- z)Q{k)) e "^-( r '- r ) 



4>{(z,T,r';Q 2 )T qq (T,b;x ¥ )T q g(r',b;x P ) . (7) 

\=L,T 

In the above, k 2 = z(1 — z)Q 2 (1 — p)/P — m 2 . Note that now, the proton is only evolved up to the rapidity ln(l/xp). 
This is because some of the energy (M x ) is carried by the dipole in order to form the diffractive final state. The 
dipole is evolved up to a rapidity ln(l//3) and the proton up to the rapidity \n{P/x) = ln(l/:rp). The high-energy limit 
in this case is ccp^C 1. 

To write formula Q, we have neglected possible final states containing gluons. This is justified because these are 
suppressed by extra powers of a s . However, if P becomes too small, the dipole evolves to higher rapidities and emits 
soft gluons. Large logarithms a s ln(l//3) coming from the emission of those soft gluons arise, and multiple gluons 
emissions should be resummed to complete formula (J2J • These multiple gluon emissions in the limit P <C 1 can also be 
accounted for in the dipole picture 0, provided one uses the large— N c limit. Indeed in this limit, a dipole emitting a 
soft gluon is equivalent to a dipole splitting into two dipoles. To illustrate this, the contribution of the qqg final state 
reads (see also 

-y* p^{X=qqg)p 
a<T diff 
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(z,|r|;Q 2 ) / d 2 z 
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T„ ( „ 2) (z,r-z,b; 



xp)-T q ,j(r,b;xp) 



(8) 



In the above, the function T qq - (z, r— z, b; xp) is the scattering amplitude for two dipoles of sizes z and r— z at impact 
parameters b— (r — z)/2 and b — z/2 respectively. These come from the splitting of the dipole of size r at impact 
parameter b. Moreover the overlap function is $^ 7 . It is so because in the leading ln(l//3) approximation, the final 
state mass Mx is fixed only by the soft gluon longitudinal momentum, and therefore transverse sizes are the same in 
the amplitude and the complex conjugate amplitude 1 . 

The diffractive vector-meson production cross-section &yM~* Vp - In this case, after scattering off the proton, 
the dipole recombines into a particular final state X = V, a vector meson whose mass we shall denote My. To describe 
this process, we need to introduce the wavefunction (z, r; My) which describes the splitting of the vector meson 
with polarization A into the dipole. The overlap function v which enters in the computation of the vector-meson 
production amplitude is then given by 



V>{' Q/3 (z,r;< 



2 )- 



(9) 



l' v \zM;Q 2 ,M 2 ) = £[^(z,r;M 2 )" 

fa/3 

Because the final state has been explicitely projected into a vector meson, the two dipoles in the amplitude and 
the complex conjugate amplitude are not connected and the cross-section (for fixed impact parameter) is simply the 
square of the amplitude: 

2 

■^ V (z,\v\ - - 



1 I'p-'Vp 
aa VM 



d 2 b 



(x v ,Q z ) = 



J d2r L dz ^, 



• Q ,M V ) T qq (r,b;x P ) 



\=L,T 



(10) 



The functions V depend on the meson wavefunctions <£>{'"^, and several different models exist in the literature 0, 
0,^3- As emphasized later, we shall only use model-independent features of V . Formula (|10|l can be used to 



compute the Deeply- Virtual-Compton-Scatering (DVCS) cross-section a 



DVCS 



provided one uses the overlap function 



< 7 (z, |r|; Q 2 ) = [^T a0 & r; 0)] * ^ a(3 (z, r; Q 2 ) 



(11) 



fa/3 

between the incoming virtual photon and the outgoing transversely-polarized real photon. 



Let us also mention another approach to compute the qqg contribution in the dipole picture, which resums logarithms of Q 2 instead of 
logarithms of 1//3. In this case, it is an effective gluonic dipole which scatters off the proton Il4l . 



III. SATURATION, GEOMETRIC SCALING AND ITS CONSEQUENCES 



Using the dipole picture of deep inelastic scattering, we have expressed the total ©, diffractive (0) and exclusive 
vector-meson production H10|) cross-sections in the high-energy limit in terms of a single object: the dipole scattering 
amplitude off the proton T qq (r, b; x). It is mainly a non-perturbative quantity, but its evolution towards small values 
of x (or high energy) is computable from perturbative QCD. Evolution equations have been established in the leading 
ln(l/a;) approximation [18i and, at least for central impact parameter, one has learned a lot about the growth of 
the dipole amplitude and the transition from the leading- twist regime T qq <C 1 towards and into the saturation regime 

T « =L " mm 

The most important prediction was probably the geometric scaling regime |2|,|20J: at small values of x, instead of 
being a function of a priori the two variables r and x, T qq (r,h ~ 0]x) is actually a function of the single variable 
r Q 2 (x) up to inverse dipole sizes significantly larger than the saturation scale Q s (x). In formulae, one can write 

T qq (r,b;x) = S(b)T(r 2 Q 2 s (x)) (12) 

where we have introduced the impact parameter profile <S*(b). Typically, S(b) = e~ h ^ R p where R p is the transverse 
radius of the proton. When performing the b integration in formulae @, or (|10|l . this contributes only to the 
normalization via a constant factor J d 2 b 5(b) = S p , caracterising the transverse area of the proton. If r 2 Q 2 > 1 then 
T =1 and the scaling is obvious. We insist that the scaling property (|12|) is a non-trivial prediction for r 2 Q 2 -C 1, 
when T is still much smaller that 1. Of course the geometric scaling window has a limited extension: at very small 
dipole sizes, deep into the leading-twist regime, the scaling breaks down. Universal scaling violations |2f)j | due to x not 
being small enough have also been derived. Recently, a new type of scaling violations has been predicted H,Ej, this 
one eventually arising when x becomes even smaller, transforming the geometric scaling regime into an intermediate 
energy regime. 

In this work, we shall consider the case of exact scaling Ijl2|l • As already mentioned, the resulti ng p rediction for the 
DIS total cross-section is in very good agreement with experimental measurements Q (see also [22|). Our goal is to 
further test the geometric scaling regime by considering its prediction for diffractive observables. But, as a reminder, 
let us start with the total cross-section. Neglecting quark masses, one can rewrite the cross-section © as 



AT poo pX 

(x, Q 2 ) = 2S P ^^ T e ) / fdf \ dz \fT(z)K 2 (^7{T^7)f) + f L (z)K 2 (^7{T^7)f)} 
11 f Jo Jo L ' 



Q 



2 



(13) 

where we have introcuced the functions /t(-z) = (z 2 + (1 — z) 2 )z(l — z) and =4z 2 (1 — z) 2 and rescaled the size 

variable |r| to the dimensionless variable f=Q|r|. We obtain the geometric scaling of the total cross-section at small 
x : 

atr X {x, Q 2 ) = v£r X (r) , r = Q 2 /Q 2 s (x) . (14) 

This has been seen confirmed by experimental data with Q s (x) given by 

-A/2 

X 



{x)=Q [ — \ , Q = lGeV (15) 



and the parameters A = 0.288 and xq = 3.04 10 4 taken from In order to illustrate it, Fig.l is an update of the 
original plot which shows the cross-section (7j f^ X as a function of r with the latest data of the different experiments 
which provide measurements at x < 0.01 : the HI ZEUS jH], E665 [H and NMC |^ collaborations. Except 
for one E665 point, the data do lie on the same curve. This is even true at low values of Q 2 , for which one could 
have exp ected scaling violations do to the charm quark mass, but one can see on Fig.l that these are not sizable (see 
also H3). 

Let us now consider the diffractive cross-section {JJ. It can be rewritten 

«" (P,x P) Q 2 ) = S p ^p^J2 e2 f dzz(l-z) fx{z)ll{z,^Q 2 s {x r )/Q 2 ) (16) 

f n P f JO \=L,T 

with the following integral 

I T , L (z,0,Q 2 s /Q 2 ) = J^fdfKr^^r^f^oWzil-zXl-^/^T^ f 2 ^j (17) 
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FIG. 1: The total cross section a 



as a function of r for x < 0.01. The data are the most recent by the HI, ZEUS, E665 



and NMC collaborations. Only statistical uncertainties are shown. 



where It contains K± and J\ bessel functions and It, contains Kq and Jo- So, another signature of the saturation 
regime of QCD should then be the geometric scaling of the diffractive cross-section at fixed /3 and small xp : 

— ^ (P,xr,Q 2 ) = — ^ (f3,r d ) , r d = Q 2 /Ql(x P ) . (18) 

This will be discussed in the next Section. Note that for small values of 0, when using a more complete formulation 
of the diffractive cross-section, the prediction (|18J) remains unchanged. For instance, in the geometric scaling regime, 

(2) 

the behavior of T^' is 

T£\*,T-z,b;x) = S(b)f(z 2 Q 2 s (x),(r-z) 2 Q 2 s (x)) (19) 

and therefore when using formula (jHJ), the prediction (|18J) remains true and it also the case in other approaches |sLll4j. 

Let us finally discuss vector-meson production. In this case, the problem is not as simple because of the extra scale 
My and of the model-dependent meson wavefunctions. However, it is possible to take advantage of a quite general 
feature rather independent of the particular model for the meson wavefunction: longitudinally-polarized vector mesons 
are predominant and the overlap function follows the approximate scaling law (see for instance |28| ) 



2nr 2 [ dz<S>l v (z,r;Q 2 ,M 2 )~g(T 2 (Q 2 +M 2 )) 
Jo 



(20) 



where the function g is sharply picked around 1. As a consequence of this, in the geometric scaling regime the 
vector-meson production cross-section can be rewritten 



cr VA/~ >yP (2*,Q 2 ) = S p 



DC- 



-q(f 2 )T QsiXr) f 2 
f 9[ ' \Q 2 + M? 



and the following behavior can be predicted: 



(21) 



4*r Vp (^ Q 2 ) - 4m* VP {tv) , r v = (Q 2 + M 2 )/Q 2 (xp) . (22) 

Note that, for the DVCS cross-section, the prediction is &'dv~cs P ( x i Q 2 ) = <t 'dv~cs P ( t ) without relying on H2U|) . 



IV. GEOMETRIC SCALING AND DIFFRACTIVE OBSERVABLES 



We shall now test the predictions 118(1 et (|22l) . The HI and ZEUS experiments at HERA have measured the 
diffractive cross section for the process ep — >eXY, selecting events with a lar ge r apidity gap between the systems X 
and Y in case of HI 2i|, and using the so-called M^-method in case of ZEUS [3(1]. Y represents the scattered proton, 
either intact of in a low-mass excited state, with My < 1.6 GeV (HI) or My < 2.3 GeV (ZEUS). The cut on t, the 
squared momentum transfer at the proton vertex, is |i| < 1 GeV 2 for both experiments. These data are presented in 

terms of the ^-integrated reduced cross section (xp, x, Q 2 ) or the diffractive structure function F^ 3 \xp,x,Q 2 ) 
obtained from the relations 



such that oy >( ' 3 ' ) —F^^ is a very good approximation, except at large y — W 2 /s with s the total energy in the e — p 
collision. HI and ZEUS measurements are realised with different My cuts so the two experiments do not measure 
exaclty the same cross-section: the proton-dissociative events are reduced in the range of the HI data set. However the 
difference is a known constant factor: ZEUS data points can be converted to the same My range as HI by multiplying 
ZEUS values by the factor 0.85 |2^,|23|. There exist also data from ZEUS for which the proton has been detected 
in the final state [3l| that we include in the following analysis. These data correspond to our definition of diffractive 
events given in Section II, as the proton truelly escapes the collision intact. But again, the obtained cross-section 
differs by only a constant factor from the one measured without tagging the final-state proton. To be comparable 
with the HI data j2£j, we multiply the data in by the factor 1.23. 

In order to test the geometric scaling properties exhibited above, we first express daJ i j^T Xp / d(3 in terms of the 
diffractive structure function: 

P = ^* P if (3) • (24) 

In Fig. 2, we present the measurements of the HI |2^| and ZEUS j^ll El collaborations for (3 da 1 di p J* Xp / d(3 as a 
function of Td = Q 2 /Q 2 {xp) at six fixed value of (3 : 0.04, 0.1, 0.2, 0.4, 0.65 and 0.90. For each of them, we include 
all data points for Q 2 values in the range [5; 90] GeV 2 and for ip < 0.01. For the ZEUS data points, the bin-center 
values in (3 given in |3JJ,|3l| are not exactly the ones quoted on Fig. 2. To be able to compare HI and ZEUS data sets, 
we have extrapolated ZEUS data points to the closest HI bin center in j3. The correction is obtained from a BEKW 
fit [s3] on the ZEUS data sets |3fi l3lj] . We have used the saturation scale 1(15(1 . It is clear on Fig. 2 that the HERA 
experimental measurements of the diffractive cross-section in DIS are compatible with the geometric scaling property 
predicted by formula 1(18(1 . as for each /3 bin, the different points form a line. 

Let us finally confront the prediction ((22(1 with DVCS and vector-mes on p roduction data from HERA. On Fig. 3, 
the available measurements from HI and ZEUS for DVCS p [23, 4> [34j and J /ip [13 exclusive productions are 
displayed. We represent the total cross-sections (meaning i-integrated) as a function of Ty (ry =t for DVCS) with the 
saturation scale ((15(1 . Again, for each vector meson, the data lie on a single curve, confirming the geometric scaling 
prediction. 



V. CONCLUSIONS 



The dipole picture of deep inelastic scattering is a theoretical framework provided by perturbative QCD (Q 2 ^> 
Xq CD ) in the high-energy limit (a;^;l). It allows to express total ©, diffractive Q) and exclusive (fj"U|) cross-sections 
in terms of a single object: a dipole scattering amplitude off the proton. One of the main features of the saturation 
regime of QCD is the scaling law (112(1 of this dipole amplitude. The resulting consequence, the geometric scaling 1(14(1 
of the total cross-section in DIS, is well-known and geometric scaling has been found in the data five years ago 
As shown in Fig.l, the present data confirm it. 

In this letter, we have exhibited the consequences of the scaling ((12(1 for diffractive observables in DIS, namely for 
the diffractive cross-section (J2J and the vector- meson production (and DVCS) cross-section ((10(1 . We have shown that 
further manifestations of the saturation scale should appear in the diffractive data in the form of the scaling laws 
((18(1 and 1(22(1 with the same saturation scale Q s as in the inclusive case ((14(1 . We have analysed the present data and 
shown them in Fig. 2 and Fig. 3. These confirm the expected behaviors and suggest that all three scaling displayed in 
Fig.l, Fig. 2 and Fig. 3 are indeed manifestations of the saturation regime of QCD. 
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FIG. 2: The diffractive cross section (3 dcrJ i F7' xp /d/3 from HI and ZEUS measurements, as a function of r<j in bins of (3 for Q 2 
values in the range [5; 90] GeV 2 and for xp<0.01. Only statistical uncertainties are shown. 



For the saturation scale Q s , we have taken over the same values of the parameters A and xq as in 4]- We did not 
try to vary them to obtain a better scaling. In the context of the exact scaling law (|12fl . a definite value for those 
parameters does not really make sense anyway. Indeed, fits of different QCD-inspired saturation models 0, on 
<jJ ^~ >X (x, Q 2 ) have shown that the precise values of the parameter are sensitive to whether or not scaling violations 
are included (they are also sensitive to what type of scaling violations are included). In any case, the parameters 
never differ significantly from those we used here. 

Note finally that in the case of vector-meson production, we only looked at cross-sections integrated over the 
momentum transfer t. Cross-sections differential with respect to t can also be expressed in the dipole picture, with the 
t dependence of the cross-section related to the impact parameter dependence of the dipole amplitude via a Fourier 
transform (see for instance USES)- We did not consider such observables in this work, but they represent natural 
places to look for geometric scaling properties at non-zero transfer as was predicted in £3] • It would require to carry 
out measurements with broad ranges for the x and Q 2 values (i.e. a large range for r), for different fixed values ott. 
This may be an experimental challenge, but it would certainly be worthy. It would especially help our understanding 
of the impact parameter (b) dependence of the dipole amplitude, and our understanding of how it is mixed with the 
high-energy evolution. 
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FIG. 3: The p, J/ty and 4> production cross-sections &vm~* Vp and the DVCS cross-section a~byc:s P from HI and ZEUS 
measurements, as a function of ry and for £p<0.01. Only statistical uncertainties are shown. 
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